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Abstract 


By comparing the observed orbital decay of the binary pulsars PSR B1913+16 
and PSR B1534+12 to that predicted by general relativity due to 
gravitational-wave emission, we are able to bound the mass of the graviton 
to be less than 7.6 x 10 _20 eV/c 2 at 90% confidence. This is the first such 
bound to be derived from dynamic gravitational fields. It is approximately 
two orders of magnitude weaker than the static-field bound from solar system 
observations, and will improve with further observations. 


I. INTRODUCTION 

General relativity assumes that the graviton has zero rest mass. For static gravitational 
fields a nonzero graviton mass m would cause the potential to tend to the Yukawa form 
r~ 1 e~ mr , effectively cutting off gravitational interactions at distances larger than the Comp¬ 
ton wavelength m -1 of the graviton. Current experimental limits on the graviton mass are 
based on the apparent absence of such a cutoff in the solar system [fj and in galaxy and 
cluster dynamics 

The study of dynamic gravitational fields allows new and independent limits to be placed 
on the graviton mass. For example, Will H and Larson and ffiscock || have shown how 
future observations of gravitational waves may be used to bound the graviton mass via 
comparison to the dispersion formula 

= \j l ~ ^2 ■ ( L1 ) 

ffere we propose a new technique for limiting m using available data on the orbital decay 
of binary pulsars. Consider the Hulse-Taylor binary pulsar, PSR B1913+16, for which the 
observed orbital decay attributed to gravitational-wave emission agrees with the predictions 
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of general relativity to approximately 0.3%. A nonzero graviton mass would alter the en¬ 
ergy emission rateQ and destroy this agreement; for gravitational waves at twice the orbital 
frequency of PSR B1913+16, requiring 

(-J < 0.003 (1.2) 

implies m < 0( 1CT 20 ) eV/c 2 . This is comparable to the limit m < 4.4 x lCT 22 eV/c 2 obtained 
from solar system observations [jlj]. 

In this paper we develop this bound in more detail. In Section |I]we examine linearized 
general relativity with a massive graviton. In Section 0 we calculate the corrections to 
the energy emission rate of a compact, slowly moving source due to a nonzero graviton 
mass. Comparison to the observed orbital decay of PSR B1913+16 and PSR B1534+12 
will provide us with an upper limit on the graviton mass in Section m We conclude in 
Section |V] with a few brief comments. 

II. LINEARIZED GENERAL RELATIVITY WITH A MASSIVE GRAVITON 

Gravitational waves on a flat background spacetime can be described as a perturbation 
of the Minkowski metric: 


9[tv T h^ u , 1 • 

For convenience we work with the trace-reversed metric, 




( 2 . 1 ) 


( 2 . 2 ) 


and raise and lower indices using r ]^, rj^. Substituting (|2.1|) into the Einstein held equations, 
expanding in powers of h^, and keeping only terms up to first order yields 


□/V = —1671-1%,, 


(2.3) 


where the stress tensor 1% of the source matter is conserved, T ,v = 0, and we work in the 
Lorentz gauge, defined by /d'A,, = 0. 

A simple generalization of ( |2.3|) to include a phenomenological mass for the graviton is 


(□ - m 2 )^ = -1671-1%,. (2.4) 

The massive held h, iv has six degrees of freedom: hve spin-2, and one scalar It is n °t 
gauge invariant, but the Lorentz condition /d%„ = 0 still holds as a constraint. 


1 Corrections to other characteristics of the system are negligible by comparison on dimensional 
grounds: ( mr ) 2 = (m /%) 2 (v / c) 2 , ( mM ) 2 = (m / lu) 2 (v / c) 6 , where v/c = O(10 -3 ) for these binary 
systems. 
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An effective stress tensor for the gravitational waves can be obtained using Noether’s 
theorem [[!(]], and is identical in form to the usual m = 0 result 

T° W = T-( W’A - • (2.5) 

Here the brackets denote an averaging over at least one period of the gravitational wave. 
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III. ENERGY EMISSION RATE 

We now calculate the energy emission rate of compact, slowly moving sources of gravita¬ 


tional waves described by (|2.4|) - (|2.5|) . Our analysis follows that used for standard linearized 
general relativity [in 


except that we work in the frequency domain for convenience |L2 


For a periodic source with period P, we decompose the metric perturbation as 




,x e 


'lOJri t 


(3.1) 


where 


n 


2tt 

~p 


(3,2) 


with n an integer. Then (|2.4j) becomes 

(V 2 + [a; 2 - m 2 ]) h^{uj\x) = -167r T^(uj\x) , (3.3) 

where V 2 is the 3-space Laplacian. The retarded Green function for this equation is 

gik\x—x'\ 


Gr(u\x-,x') = 


(3.4) 


47r|T — x'\ ’ 

where k = sign(u;)\/co’ 2 — m 2 for |o;| > m. The retarded solution of (|3.3|) for fixed u is then 


hfj, v {ui\x) = 167T / d 3 x' Gr(u\x; x') \x') . 


(3.5) 


Using the slow-motion approximation (uia <C 1, with a the characteristic size of the source), 
taking the observation point far from the source region (r = |x| 3> a), and employing the 
conservation of the stress tensor, one obtains 


h 0 o(uj\x) 
h 0 j(u\x ) 
h jk (u\x) 


4 ,gikr 

r 

4 ,gikr 

r 

4, a ikr 

r 


„ ,—. 'f-.y rv> k —, 

M + -(-ik)D J + —(-ikyi jk 
—{—iu])Dj — —(—ik)(—iw)T jk , 
iu) 2 I jk , 


(3.6) 
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where M , Dj , are respectively the Fourier coefficients of the mass, dipole moment, and 

quadrupole moment of the source. Only the quadrupole terms contain nonzero-frequency 
components and contribute to the radiation. 

Substituting ( |3.6|) into the effective stress tensor (|2.5| ) for the gravitational waves and 
integrating the outward flux over a sphere centered on the source gives the rate of energy 
emission. One finds 


L~ dE — L 

L = ~n~ Lcr 


+ £ 


00 2 4 

m oj„ 


n= 1 


Ijk (^n ) Ijk (Ai) tr I (u n ) 


+ O (m 4 ) , (3.7a) 


where 


l gk = 




71=1 


7/j7 C (cu n )/.^.(tu ri ) tr /(ttJyj) 

' 15 


(3.7b) 


is the usual general-relativistic expression for the radiated power, tr / is the trace of Ijk, and 
we sum over repeated indices. Equation Q3.7a|) gives the corrections to the radiated power 
due to a small nonzero graviton mass; comparison to the observed orbital decay in binary 
pulsars PSR B1913+16 and PSR B1534+12 will provide us with a bound on m. 


IV. BINARY PULSARS 

The formula (|3.7| ) for the energy-loss rate of a gravitational-wave source when the gravi¬ 
ton is massive is easily applied to the orbital decay of binary systems to put a limit on 
the graviton mass. Consider PSR B1913+16, for which the orbital decay rate is slightly in 
excess of the predictions of general relativity || . Denote by p, the measured orbital period 
of the binary system, Pb the measured orbital period derivative ascribed to gravitational 
radiation, and Pgr the instantaneous period derivative expected from general relativity. For 
a slowly decaying Keplerian binary, the instantaneous period derivative is proportional to 
the energy-loss rate; hence, 


A = 


Ph-P 


GR 


L — L 


GR 


P 


GR 


J GR 


(4.1) 


where L is the gravitational-wave luminosity inferred from Pb, and Lgr is the energy-loss 
rate expected from general relativity. This fractional discrepancy A has been measured for 
PSR B1913+16 and PSR B1534+12 (see [|6|.ii7] and Table [I]). 

Now suppose that A is due at least in part to a nonvanishing graviton mass. Combining 
(fh7[) and (|Q|) implies 


24 

m 2 < — P(e) 


2irh \ Pb — P< 


c 2 P, 


GR 


P 


(4.2) 


GR 


where P(e) is a function of the eccentricity, 
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(4.3) 


E 


n 


F{e) = 


1 n— 1 


3 1jk^^ri) tl I(u)ri) 


12 


E 

n=l 


Ijk (^n) 11 ' / ( UJ n ) 


1 + 73 2 37 4 

A 24 c 2! 96 

(1 — e 2 ) 3 


as can be shown using the techniques of |L3|; see Figure |T]. 

Equations ( |4.2p , ( |4.3| ) show that a nonzero graviton mass increases the energy emission 
and decay rate of Keplerian binaries, as one would expect from adding extra degrees of 
freedom to the gravitational field. The strongest bounds arise from binaries with small 
eccentricity and large period, as these systems produce the bulk of their radiation at low 
frequencies [|R|, which are the most sensitive to a graviton mass, as in equation ( jl.l|) . 

In using ( |4.2| ) to place an upper limit on the graviton mass, we should take into account 
the experimental uncertainties in the fractional discrepancy A, which are typically of the 
same order as A. We assume the measured discrepancy A to be normally distributed 
about its unknown actual value (given by the equality in (|4.2|) with unknown m 2 ), and 
with standard deviation as given in Table In our model we relate the discrepancy to the 


squared graviton mass, which must be non-negative. Referring to Jl4], Table X], which lists 
the 90% unified upper limit/confidence intervals for the non-negative mean of a univariate 
normal distribution based on a measured sample from the distribution, we calculate the 
90% upper limit on the graviton mass, which is given in the final row of Table |. These 
two observations of m 2 may also be combined into a single upper bound on the graviton 
mass by averaging the individual m 2 bounds with weight according to their variances. Again 
referring to Table | and ||T4| , Table X], the corresponding limit on the graviton mass from 
the combined observations of PSR B1913+16 and PSR B1534+12 is found to be 


ttt. 90 % < 7.6 x 10 20 eV/c 


(4.4) 


V. DISCUSSION 

Table | gives the relevant parameters and the corresponding graviton mass bounds for 
the two binary pulsars whose gravitational-wave induced orbital decay has been measured, 
PSR B1913+16 and PSR B1534+12 [[].[7|. These bounds are about two orders of magnitude 
weaker than the Yukawa limit obtained from solar-system observations, me 2 < 4.4 x 10~ 22 eV 
II, and several orders of magnitude weaker than that provided by observations of galactic 
clusters, me 2 < 2 x lCT 29 eV P||, though these galactic cluster bounds may be less robust, 
owing to their reliance on assumptions about the dark matter content of the clusters. In 
contrast, the bound obtained here is very straightforward. Our chief assumption is the form 
of the effective mass term for the graviton, which, while not unique, is natural. Furthermore, 
any other mass term would be expected from dimensional arguments to yield similar results. 

Our other major assumption is that only unbiased measurement errors enter into the 
determination of the intrinsic binary period decay rate P\j. The determination of P f, requires 
an accurate distance measurement to the binary system, however, which can be difficult 
to make. The large uncertainty in the discrepancy A associated with PSR B1534+12 may 
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well be due to an underestimate of the distance to this binary system 0, in which case the 
bound on m 2 would be even tighter. 

The bound described here arises from the properties of dynamical relativity, making 
it conceptually independent of either the solar system or galactic cluster bounds on the 
graviton mass, which are based on the Yukawa form of the static field in a massive theory. 
Furthermore, the bounds from any given pulsar system will improve as observations increase 
the accuracy of the measured fractional discrepancy in the period derivative. 
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TABLES 


TABLE I. Orbital parameters and corresponding graviton mass bound from the two binary pul¬ 
sar systems whose gravitational wave induced orbital decay has been measured. Pulsar parameters 
are taken from [||,j7]]. One-signra uncertainties are quoted for A. 



PSR B1913+16 

PSR B1534+12 

Period 

27907 s 

36352 s 

Eccentricity 

0.61713 

0.27368 

A 

0.32% ± 0.35% 

-12.0% ± 7.8% 

Graviton mass 90% upper bound 

9.5 x 10 -20 eV/c 2 

6.4 x 10~ 20 eV/c 2 








FIGURES 



0 0.2 0.4 0.6 0.8 1 


e 

FIG. 1. Eccentricity factor F(e) (cf. eqn. [T^) versus e. 
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